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1. Motivation 

In order to determine renormalized quantities like quark masses, the chiral condensate or 
the pion decay constant on the lattice one needs renormalization factors to compare with 
results from experiments or continuum theories. Those values are typically given in the 
modified minimal subtraction (MS) renormalization scheme. 

Lattice Dirac operators satisfying the Ginsparg- Wilson (GW) equation Q implement a 
version of chiral symmetry that is closest to the continuum form, with only local O (a) 
violations. Exact GW-operators have the benefit of relations that simplify the renormal- 
ization procedure 0] . For Dirac operators that fulfill the GW-equation only approximately 
computing renormalization factors also provides a check how well chiral symmetry is real- 
ized in this formulation. 

Let us consider local, flavor non-singlet bilinear quark-field operators 

O r = uTd, (1.1) 

where T denotes a matrix from the Clifford algebra in the chiral representation and u 
and d denote the quark fields for the light quarks. We organize the 16 matrices into the 
scalar, vector, tensor, axial vector and pseudoscalar sectors with S ~ t,V ~ 7/njT ~ 
\ [7/LtjTf] ~ 7//Y5, P ~ 75 according to their properties under Lorenz transformations. 
For a chiral Dirac operator we have Zp = Z$ and Za = Zy . For conserved covariant 
currents we have Za = 1 due to Ward identities. 
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Although only the overlap operator ||, Q satisfies the GW-equation exactly, there is a 
number of formulations that have good chiral properties by satisfying the GW-constraint 
approximately (fixed point operator ||, Dq\ §) or in some limit (domain wall fermions 
(7[ ^]). Within the BGR-collaboration fixed point and CI operators have been studied in 
both, the quenched approximation || and the CI operator was analyzed in full QCD in 

The Dei has been introduced in B, |l3| as an approximate solution to the GW-equation, 
where the lattice derivative operator does not only include nearest neighbor interactions, 
but also more remote connection, each of which carries a Dirac structure. 



Renormalization factors for the quenched case have been computed in [14| and this work 



follows the procedure therein closely. Non-perturbative renormalization in the quenched 



approximation has also been performed for the Wilson operator [15, 1£, 17, [Tj|], staggered 



fermions [19|, domain-wall fermions [£(J] and the overlap operator [21, [2^, |23|]. Calculations 
with dynamic fermions were performed non-pertubatively for the overlap operator in [24]. 
For dynamical clover fermions perturbative results were presented in [^, 26], while [27] 
covers perturbative and non-perturbative approaches. 

In Sect. ^| the method to determine renormalization constants [28] is reviewed. Technical- 
ities like the parameters used for generating gauge configurations and data analysis are 



dealt with in Sect. y. Results in different schemes are discussed in Sect. [4.1| , Sect. and 
Sect. and we summarize our results and conclude in Sect. ||. 

2. Method 

We want to compute non-perturbative renormalization constants on the lattice using a 
scheme applicable in both, lattice simulations and perturbative continuum calculations. 



To this end we use the regularization independent scheme (RI) proposed in [28]. Results in 
the Rl-scheme can also be converted to the more conventional MS-scheme in a perturbative 
way. In the Rl-scheme we match expectation values of bilinear quark operators between 
quark fields at a specific momentum value p 2 = fj, 2 with the renormalization scale /i, 

(p\Or\ P )\ p 2 =fl , (2.1) 

to corresponding tree-level matrix elements {p\Or\p) - Continuum perturbation theory 
can only be applied for a renormalization scale much larger than the QCD scale parameter 
Aqcd- Discretization effects can be neglected for fi much smaller than the lattice cut-off 
1/a so the renormalization procedure is expected to work in a window 

A^cd « f « I/" 2 • (2-2) 

Expression fl2.1| ) is gauge-variant so we need to fix the gauge in order to compare results 
obtained in a lattice formulation to continuum results. Landau gauge can be implemented 



in both approaches, but has the problem of Gribov ambiguities. In |29| , p0| , 31, 14] this 
has been addressed and no significant effects were detected, therefore no explicit check was 
performed in this work. 
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We use the method from [28] with the modifications from |32| to compute renormalization 
factors. Multiplying (2.1) with the inverse of the tree-level matrix element (p|Or|p)o 1 gi yes 



us the renormalization condition 



Z T — Tr 
12 



-i 



The matrix element 



(p\O r \p) (p\O r \p) 



(p\O r \p) = ^A r (p) 



(2.3) 



(2.4) 



with the quark field renormalization factor Z q is proportional to the amputated Green 
function 

Ar(p) = S- 1 (p)G r (p)S- 1 (p). (2.5) 

The inverse quark propagator is denoted S~ 1 (p). The Green function Gr(p) is computed 
by taking the expectation value of the local operator Or (z) and Fourier transforming it. 
This reads 



Gr(p)a,t 



v ^ 



e ip{x ~ y) ( u r 



(x)J20r(z)dp(y)) , 



(2.6) 



where the indices a, f3 run over color and Dirac indices and V denotes the lattice volume. 
We use equal masses for u and d quarks, hence the quark propagator is 

S a ,fs(x,y) = (u a (x)up(y)) = (d a (x)dp(y)) . (2.7) 

Using the modification (3^] we compute quark propagators with momentum sources 

S n (x\p) = J2^ PV Sn(x,y) (2.8) 



for each gauge configuration. Inserting the local operator defined in Eq. ( |l.l| ) into Eq. ( [2.6] ) 
we can rewrite Gr(p) in terms of the quark propagator S(x\p). Using 75-hermiticity, we 
rewrite the Green function 



-ip(x-y) 



u a (x)up(z)Td a (z)dp(y)) 



x,y,z 
N 



^T7 E E l5Sn{z\p) j ^S n (z\p) , 



(2.9) 



VN 



n=l z 



where N is the number of gauge configurations we average over. The quark propagator in 
momentum space is obtained by transforming Eq. (p.8|) fully to momentum space 



-ipx 



S n (x\ 



(2.10) 



71 = 1 X 



We use momentum sources [32] with momenta listed in Table ^ and Table |3| in order to 
compute S n (y\p) by solving 



^D cl (x,y)S n {y\p) = e 



ipx 



(2.11) 
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for the chirally improved Dirac operator Dqi- This implies we have to solve Eq. ( [2.11 ) for 
every momentum vector, but improves the signal significantly. In momentum space this 
Dirac operator reading 

D cl (p) = R{p) + i 1 ^(p) (2.12) 



in the free massless case is constructed using the functions 0, 13] 



a ^(p) = iPti + O {ap) 2 and R(p) = O (ap) 2 . 



(2.13) 



The quark field renormalization is obtained in the so-called Rl'-scheme, where we can apply 
Eq. ( [Hp and find 



Zf\p) = ±-TrS-\p)S {p) 



—Tr /V 1 ^ R (P)- ia »(PH 
12 V yP >R(p) 2 + a,{p)a,(p) 



(2.14) 



Renormalization factors for quark bilinears can then be computed from the renormalization 
condition Eq. (|2.3|) together with the definitions Eq. (|2.4j) and Eq. (|2.5|) for the Rl'-scheme 



Zf l \p) 



12 



Tr 



(p\O r \p) {plOrlp)- 1 
\2Zf\p) 



Tr [S- 1 (p)G T (p)S- l (p)T 



(2.15) 



where we use (p | Or | p) = T. For the vector, axial vector and tensor sectors we average 
over the components before taking the trace. 

A modified version of the method described above was presented in |?3| and also in |34j , 
which will be referred to as the "reduced" method. The contribution that comes with the 
unit matrix of the quark propagator introduces a pure cut-off effect and should therefore 
be eliminated from the renormalization prescription. A way to achieve this is to define a 
reduced quark propagator 



- ■%,) - [^TrS(p)} i 



(2.16) 



as well reduced free quark propagator 

So(p) = S (p) - 



TrS (p) 1 



(2.17) 



Consequently a reduced quark field renormalization factor 

z T (p) = -^ Tr ~s 1 (p) s o{p) 



(2.18) 
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can be computed. Analogously Eq. ( |2.15| ) can be redefined in these terms 



Zf\p) 



I2zf{p) 



Tr 



S \p)G T {p)S \p)T 



-i 



(2.19) 



with 



1 N 

Grip) ~ yjj Yl Yl l5S„(z\p)^ 5 rS n (z 



{ z\p) . (2.20) 

n=l z 

We state the final numbers for the renormalization factors in both, the non-reduced and 
the reduced scheme. 



3. Technicalities 
3.1 Parameters 

Gauge configurations were generated with two dynamic quarks using the Liischer-Weisz 



gauge action [35] and stout smearing [36], which are presented in [10, 11, 12]. Prior results 
in the quenched approximation utilized hypercubic smearing instead of stout smearing, 
hence a comparison of renormalization factors for the quenched and the dynamic case are 
not applicable. The gauge couplings and sea-quark masses are listed in Table [j]. For the 
12 3 x 24 setups we compute 10, for the 16 3 x 32 setups 5 momentum propagators for each 
momentum vector. 



run 


L 3 x T 


/3 


a m q 


omAwi 


m Awl [MeV] 


a[fm] 


alGeV" 1 ] 


#cf. 


a 


12 3 


x 24 


4.7 


-0.05 


0.023(2) 


30.9(8) 


0.147(18) 


0.75(9) 


10 


b 


12 3 


x 24 


5.2 


0.02 


0.025(1) 


42.1(2) 


0.115(6) 


0.58(3) 


10 


c 


12 3 


x 24 


5.2 


0.03 


0.037(1) 


58.1(2) 


0.125(6) 


0.63(3) 


10 


d 


12 3 


x 24 


5.3 


0.04 


0.037(2) 


60.8(2) 


0.120(4) 


0.61(2) 


10 


e 


12 3 


x 24 


5.3 


0.05 


0.050(19) 


76.4(2) 


0.129(1) 


0.654(5) 


10 


f 


16 3 


x 32 


4.65 


-0.06 


0.02664(15) 


35.05(19) 


0.150(1) 


0.760(5) 


5 


g 


16 3 


x 32 


4.70 


-0.05 


0.03277(31) 


43.11(41) 


0.150(2) 


0.76(1) 


5 


h 


16 3 


x 32 


4.58 


-0.077 


0.01097(26) 


15.03(36) 


0.144(2) 


0.73(1) 


5 


i 


16 3 


x 32 


4.65 


-0.07 


0.00861(36) 


12.09(52) 


0.140(1) 


0.711(5) 


5 



Table 1: Parameters of the different runs with AWI mass and lattice spacing measurements thereon. 
The lattice spacing and AWI mass values are taken from jl2], |l0| . For run i the values are preliminary 
and not yet published. Runs a — e denote runs on the smaller 12 3 x 24 lattices, while runs / — i 
denote the runs on the 16 3 x 32 lattices. Note that the lattice spacing a for run a is closer to runs 
with big lattice size than to the remaining four 12 3 x 24 lattices. 



3.2 Interpolation, pole contributions and chiral limit 



In the quenched case [14] we performed a chiral extrapolation of the renormalization factors 
using the different valance quark masses for each lattice setup and then stated these values 
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1 


U.ylDo 


l.zoUU 


1.0 f ZO 


1 A AdZ 
1.440D 


l.oUuo 


1.4U1D 


1 


1 


1 


1 


9883 


1.3266 


1.6958 


1.5601 


1.6251 


1.5117 


2 


1 


1 


1 


1.1184 


1.5013 


1.9191 


1.7655 


1.8391 


1.7108 


3 


1 


1 


1 


1.2892 


1.7306 


2.2121 


2.0352 


2.1200 


1.9721 


2 


1 


1 


2 


1.4399 


1.9329 


2.4707 


2.2730 


2.3678 


2.2026 


3 


1 


1 


2 


1.5762 


2.1159 


2.7047 


2.4883 


2.5920 


2.4111 


4 


2 


2 


2 


2.1628 


2.9033 


3.7112 


3.4143 


3.5565 


3.3084 



Table 2: The available momentum combinations on the 12 3 x 24 lattice. Parameters of the runs 
are listed in Table [j]. 



api 


ap2 


OP3 


aPA 


ap 


f 

p[GeV] 


g 

p[GeV] 


h 

p[GeV] 


i 

p[GeV] 


2 


1 


1 


1 


0.8388 


1.1011 


1.1035 


1.1494 


1.1802 


3 


1 


1 


1 


0.9669 


1.2693 


1.2720 


1.3250 


1.3604 


2 


1 


1 


2 


1.0799 






1.4798 




4 


1 


1 


1 


1.1151 






1.5280 




3 


1 


1 


2 


1.1822 


1.5519 


1.5552 


1.6200 


1.6633 


4 


1 


1 


2 


1.3061 


1.7146 




1.7898 




3 


1 


2 


2 


1.3639 


1.7904 


1.7942 


1.8690 


1.9189 


3 


2 


2 


2 


1.5241 


2.0007 


2.0049 


2.0885 


2.1443 


4 


2 


2 


2 


1.6221 


2.1294 


2.1339 


2.2228 


2.2822 


6 


1 


2 


2 


1.7369 


2.2801 


2.2849 


2.3801 


2.4437 


5 


1 


2 


3 


1.8235 


2.3938 


2.3989 


2.4988 




5 


2 


2 


3 


1.9462 


2.5549 


2.5603 


2.6670 


2.7383 



Table 3: The available momentum combinations on the 16 3 x 32 lattice. The p values are only 
listed if we computed the quark propagator for this momentum vector. Parameters of the runs are 
listed in Table [j]. 



as our final numbers for Zr- In a full QCD calculation we only have one physical quark 
mass per lattice setup, so a chiral limit cannot be performed in a straightforward way. 
In this work we present the renormalization factors for each lattice setup and perform 
an extrapolation using different valence and sea quark masses, i. e. a partially quenched 
approximation, only for discussing the continuum limit in Sect. iA and for extracting the 
pion pole contributions. The pseudoscalar density couples to the Goldstone boson channel, 
so we expect O (1/m) contributions to the pseudoscalar propagator in the chiral limit. As 
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run 


a rriq 


am q 


a 


—0.05 


-0.025 


0.0 


0.03 


0.06 


0.10 










b 


0.02 


0.065 


0.11 


0.155 


0.20 












c 


0.03 


U.U I ZO 


n iir 


U.lD 1 


n on 
U.zU 














U.U1 


0.08 


0.12 


0.16 


0.20 












e 


0.05 


0.875 


0.125 


0.1625 


0.20 












f 


-0.06 


0.0 


0.06 
















g 


-0.05 


0.0 


0.05 
















h 


-0.077 


-0.067 


-0.057 


-0.047 


-0.037 


-0.027 


-0.017 


0.003 


0.026 


0.063 


i 


-0.07 


-0.06 


-0.05 


-0.04 


-0.03 


-0.02 


-0.01 


0.0 


0.02 


0.04 



Table 4: Valence quark masses used in the configurations in Table [l| The second column denotes 
the case where sea quark mass and valence quark mass are equal, i. e. full dynamic calculations, while 
the columns three to eleven mark partially quenched situations used for the pion pole extractions 



in Sect. 3.2 



discussed in 37, 39, 40 1 the inverse renormalization factor can be expanded as 



— = — +B + Cm + (m 2 ) 
Zp m ' 



(3.1) 



and a "subtracted" renormalization factor can be defined by subtracting the pole term 



1 



1 



Z% uh Z P 



A 



in 



(3.2) 



The operator product expansion ensures that pole contributions are suppressed for large 
values of [i. For runs a — e,h and i we have sufficiently many different valence quark masses 
and we use the lowest five to perform a fit to Eq. fl3.1|) . For runs / and g we only have three 
valance quark masses, so here we only solve Eq. ( |3.1[ ) for the coefficient. The valence quark 
masses are collected in Table ||. Fig. [I] shows the original and subtracted renormalization 
factors Z P and Z s P uh of runs h and i for selected momentum vectors plotted against the 
AWI mass including the corresponding fit. 

In Fig. § the coefficient A of such a fit is displayed for runs a — e, h and i. The coefficient 
A is expected to be proportional to {qq) /p 2 , where — (qq) = S is the chiral condensate 
(cf. In one- loop partially quenched chiral perturbation theory the condensate 

includes a logarithmic term that leads to a divergence in the chiral limit [^2| , but a possible 
additional logarithmic term cannot be disentangled from the 1/m pole term with valance 
quark masses available. 1 Coefficient A shows a rather universial behavior with a estimated 
uncertainty of 3% for (qq), which indicates a consistent pion-pole removal over the whole 
range of lattice configurations. In order to compare runs at the same physical momentum 
transfer we use a cubic square interpolating fit. The values presented at momentum transfer 
/i = 2 GeV are taken from such interpolating fits. 

Uncertainties of all quantities are calculated using the standard jackknife procedure. Note, 



1 Thanks to the referee for pointing this fact out. 
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however that we include only uncertainties of statistical nature and no systematics due to 
finite size effects or other artifacts. 

For selected momenta we have higher statistics available for run a (23 configurations) and h 
(15 configurations) from Table |]. As expected uncertainties are lower for larger ensembles, 
but in the region with > 1 GeV all mean values for the renormalization factors with 
larger sample number lie within the error band of the results with smaller sample number. 
Only in the region below 1 GeV we find larger deviations, but those values also come with 
large uncertainties. 
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Figure 1: The renormalization factors Zp and Zp plotted against the the AWI-mass for selected 
momenta of runs h and i. A fit to Eq. (3.1) is shown for Zp and the same fit with subtracted pole 
term is plotted for Z§, ub . 



4. Results 

4.1 Results in the Rl'-scheme 

In Fig. H and Fig. |4] we present the result of Z Rr plotted against the momentum for runs on 
12 3 x 24 and 16 3 x 32 lattices, respectively. Additionally an interpolating fit is displayed. 
Runs b — e seem to agree for most of the quantities, whereas run a lies closer to the runs 
on the 16 3 x 32 lattices. Therefore run a is also included in the latter plot for comparison. 
This separation might be a hint for finite size effects, because run a has a significantly 
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Figure 2: The coefficient A of the fit to Eq. ( |3.f[ ) displayed against the /z 2 . We plot the runs a 
(black circles), b (red squares), c (blue diamonds), d (green upwards triangles), e (brown downwards 
triangles) on the 12 3 x 24 lattices and the runs h (violet left pointing triangles) and i (cyan right 
pointing triangles) on the 16 3 x 32 lattices. The parameters of the different runs are presented in 
Table 0. 



larger box size then the rest of the 12 3 x 
region of runs f — i. 

A chiral operator is expected to have Za 



24 setups and also the lattice spacing lies in the 



end we also plot the ratios Zp Ub,Rl /^Ip 



Zy = const, and Zp = Z$ = const.. To this 
s . - Z S p uh jZ s and Z^ v /Zy 1 ' in Fig. g, where the 
scaling behavior is supposed to cancel out. In the interval 1.5 GeV < \i < 3.0 GeV this 
is satisfied up to 7% and 10% for Z p ub ' Rl /Zg 1 and ZY /Zy 1 = Za/Zv, respectively, on 
the larger lattices. The Dqj does not satisfying the GW-relation exactly, so we cannot 
expect the ratios to be equal to unity. In the quenched case [14] a comparable behavior 
was found for Za/Zv, while Zp ub /Zs was closer to unity. Note, that these ratios are the 
only quantities that we can directly compare to other studies, as the absolute values of the 
renormalization factors strongly depend on the formulation of the Dirac operator and the 
smearing method in use. 



4.2 Results in the MS-scheme 

The difference between the Rl'-scheme we use and the Rl-scheme lies solely in the definition 
of the quark field renormalization factor, hence 

Zf _ Zf l 



JtPRI 



7RI' 7 RP 
Zj q Zj v 



(4.1) 



The conversion factor expanded in terms of the strong coupling constant a s in Landau 



gauge reads [ 43 1 



cf RI =1 + 



_67 2^1/0^2 

6 3 J Uvr/ 
52321 607C 3 /2236 

+ { -— + -r + W 



*KW£)' + °(£) 4 . 



(4.2) 



where Nf denotes the number of dynamic quarks and C, n denotes the Riemann zeta function 
evaluated at n. 

Due to Ward identities the conversion factor for axial vector and vector operators reads 



JU'MS _ RI'RI RIMS _ JU'RI _ JU'RI : JIIMS _ -1 
C V,A — C F C V,A — c r — c q L e - C V,A — 1 ■ 



(4.3) 
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Figure 3: The renormalization factor Z m ' displayed against /i 2 for all 12 3 x 24 lattices 



For scalar and pseudoscalar sectors we have [44] 
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Figure 4: The renormalization factor Z m displayed against /j, 2 for 12 3 x 24 with j3 = 4.7 and all 
16 3 x 32 lattice setups. 



From |^5|, ^] we get the conversion for the tensors 
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I88OC3 80C 4 225890 \ 



27 



+ 



+ 



729 



-jN f + 



/_32C3_ 9542\ 2 Wo,y 



(4.5) 
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q 



10 



\x 2 [GeV 2 ] 



Figure 5: The ratios Z§, ub /Z s and Z^/Zy vs /i 2 for 12 3 x 24 with /3 = 4.7 and all 16 3 x 32 lattice 
setups. 



Finally for the quark field we use [46] 



„RIMS 



( 1622 184C 3 2597V/ ] /a.y 
\ 27 9 81 J Utt/ 



+ 



15479317 1209445C 3 1072C 4 10040C 5 



5832 972 81 

/ 225890 I88OC3 8OC4 



729 



27 



+ 



9 



N f + 



27 



2187 81 J / J V4vr/ V4vr/ 



The 3-loop expression for the coupling a s in the MS-scheme [47] reads 

«s(g 2 ) 1 _ A log(log(g 2 )) 

4tt /3 log((Z 2 ) /3 3 log(g 2 ) 2 

+ 05 1 \ 2,3 W lo g(M<Z 2 )) 2 " ft 2 log(log(9 2 )) + ~ Pi) 
f3^\og(q 2 ) 6 



(4.6) 



(4.7) 
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with q = ///Aqcd- The QCD scale (from [48] with ro = 0.5 fm) in the MS-scheme is given 

by 



A^" 2 = 0.243(24) GeV (4. 



and we use the following coefficients 



A) = 11 - \n s , (4.9a) 
ft = 102 - —N f , (4.9b) 



3 

2 ~ ~18 " J ' 54 



, m _ 2857 5033 . . , 325 - 



?2 = — ~ N f + ~ N f > ( 4 - 9c ) 



MS 149753 , Z' 1078361 6508Cs 



«"•=— +3BM<a+^^ 27* ) A > 

50065 64720A ^ 1093^ 
162 81 / y 729 i ' 



(4.9d) 



The resulting renormalization factors in the MS-scheme at fj, = 2 GeV are collected in 
Table [| and Table ||| for the original and the reduced definition, respectively. 

4.3 Results in the RGI-form 

Renormalization factors in general depend on the renormalization scale jx in a way that is 
determined by the anomalous dimension 70 of the operator in consideration 

lo(a s ) = -^lnZo(/x 2 ) = f> • (4-10) 



i=0 

After integrating Eq. ( |4.10|) we arrive at 

p\ \ ~7o/(2/3o) 

4?r 

the renormalization factor in RGI-form (Renormalization Group Invariant), where the scale 
dependence coming from the renormalization group is extracted up to a certain order in 
perturbation theory. We use the the 4-loop expression of the QCD-/3-function 




z RGI = z(f) ( 2A^4^ I 'Mil ''"I — + — U- 11.111 



? = = Eft©'". <«■*> 



with coefficients listed in Eq. (4.9). Axial vector and vector are scale independent hence 
their anomalous dimensions are 7^4 = jy = and therefore Z^y = Z~^y- For j q , 7s = 
7p = 7" 1 and 77 1 we use the 4-loop expressions from [49, 51, |45[ |. 

In Fig. ^ and Fig. ^ the renormalization factors in RGI-form are displayed against the 
momentum transfer \j? for the 12 3 x 24 and 16 3 x 32 lattices, respectively. In Fig. |8| the 
reduced renormalization factors in RGI-form are presented against the same quantity. In 
the interval 1.5 GeV < fi < 3.0 GeV we find a maximal deviation from the plateau behavior 
of 8% for Zs and 5% for the remaining renormalization factors. In the case of the reduced 
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factors Zq, Zy, Za and Zt we observe an almost linear scaling behavior in the region 
[i > 1.55 GeV or a/j, > 1.15. A linear fit of the form c(l + d(afi) 2 ) allows us to estimate 
discretization errors proportional to (afj,) 2 . For the afore mentioned quantities we find slope 
values d of 0.019, 0.042, 0.020 and 0.039. 






ji 2 [GeV 2 ] 

Figure 6: The renormalization factors in RGI-form Z RGl plotted against for the 12 3 x 24 
lattices, runs a — e in Tabic 
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4.4 Collection of Results 

In this section we collect results of renormalization factors for the original and the reduced 
method as well as for massive quarks and a chiral extrapolation in Tables |^-||. The same 
quantities are plotted in Fig. || against the lattice spacing. The difference between the 
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original and the reduced method lies in the subtraction of the component of the quark 
propagator that is proportional to the unit matrix for the latter method (cf. Eq. ( 2,16| )). 
Thereby a potential cut-off artefact is reduced according to 1 33 J- It is interesting to note that 
the renormalization factors scalar, vector, tensor and axial vector in the reduced definition 
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7 Sub 
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n (111 fJ/Q*! 

U.yiio(o J 


u.y i o(4 j 


U.oo40(o ) 


U.yuy4( i ) 


u.yio(i j 


U. ( t y(o j 


43(2) 


n 1 1 

U. 110(0 J 


U.yo ( z(o J 


u.yo(i J 


U.oyo^l J 


U.y /o(l J 


U.yoo(z J 


U.ooo(4 J 


Oo(o ) 


U.lzo(0 J 


U.yooy( r J 


u.yy < (4j 


U.oyD^l J 


U.y / z(l J 


U.yt>t>(z J 


U.o4y(z J 


61 (9\ 


190(4"l 


9739(9"l 


n 989(5") 


900(9"l 


978(1 "l 


968(1 "l 


850(3") 


76.5(6) 


0.129(1) 


0.9763(3) 


1.021(2) 


0.9060(6) 


0.9758(5) 


0.9739(4) 


0.881(2) 


35.0(3) 


0.150(1) 


0.8987(5) 


0.969(3) 


0.8257(9) 


0.8967(5) 


0.902(1) 


0.783(2) 


43.1(5) 


0.150(2) 


0.9030(3) 


0.972(4) 


0.8290(9) 


0.8998(4) 


0.9068(5) 


0.800(4) 


15.0(4) 


0.144(2) 


0.8934(7) 


0.93(1) 


0.818(2) 


0.896(2) 


0.894(2) 


0.71(1) 


12.1(5) 


0.140(1) 


0.9013(3) 


0.93(2) 


0.826(2) 


0.9042(8) 


0.901(1) 


0.707(8) 



Table 5: The interpolated rcnormalization factors in the MS-scheme Z ° with the corresponding 
AWI mass toawi & t momentum transfer fi — 2 GeV. 



m AW i[MeV ] 


a[fm] 


z q 


Z s 


Zy 


Zj< 


Z A 


ySub 


31(4) 


0.147(18) 


0.961(1) 


1.025(5) 


0.8798(9) 


0.959(1) 


0.964(2) 


0.855(3) 


43(2) 


0.115(6) 


1.005(2) 


0.97(1) 


0.928(2) 


1.014(2) 


0.995(3) 


0.912(5) 


58(3) 


0.125(6) 


1.027(2) 


1.057(4) 


0.950(2) 


1.030(2) 


1.023(2) 


0.930(2) 


61(2) 


0.120(4) 


1.031(1) 


1.047(5) 


0.953(2) 


1.035(2) 


1.024(2) 


0.929(3) 


76.5(6) 


0.129(1) 


1.0477(7) 


1.096(2) 


0.9722(9) 


1.047(1) 


1.0452(9) 


0.976(2) 


35.0(3) 


0.150(1) 


0.9504(5) 


1.025(3) 


0.8732(9) 


0.9483(6) 


0.954(1) 


0.846(2) 


43.1(5) 


0.150(2) 


0.9604(4) 


1.033(4) 


0.8817(7) 


0.9570(4) 


0.9645(6) 


0.869(4) 


15.0(4) 


0.144(2) 


0.9278(6) 


0.96(1) 


0.849(1) 


0.930(1) 


0.928(2) 


0.75(1) 


12.1(5) 


0.140(1) 


0.9327(4) 


0.96(2) 


0.855(2) 


0.9358(9) 


0.932(1) 


0.740(7) 



Table 6: The reduced, interpolated renormalization factors in the MS-scheme Z with the cor- 
responding AWI mass toawi a t momentum transfer /i = 2 GeV. 



differ from the original definition solely by a factor Z q jZ q before taking the chiral limit, 
which means that the denominator in Eq. ( |2.19| ) is not influenced by the redefinition of 
the quark propagator. For the pseudoscalar renormalization factor we observe a maximum 
deviation of 5% from this behavior. For chirally extrapolated values this strict factorization 
in no longer observed, but still holds approximately, with deviations in the sub percent 
range for vector, tensor and axial vector, and below 3 % and 5 % for the scalar and 
pseudoscalar, respectively. 

5. Summary and Conclusion 

Renormalization factors are essential to relate computations of renormalization scheme de- 
pendent quantities like the pion decay constant and quark mass on the lattice with results 
from continuum calculations or measurements. We present values of the renormalization 
factors for quark bilinears in the regularization independent scheme and the modified min- 
imal subtraction scheme using conventions of |28|, ^] and [^] . For the pseudoscalar factor 
the pion pole is subtracted and a chiral extrapolation was performed. 
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r tv yr T 7 1 

m A wi[MeV J 


a[tmj 


z q 


Zs 


ry 

Zy 


ry 

ZjT 


Za 


rySllb 


31(4) 


0.147(18) 


0.9007(3) 


0.948(4) 


0.8220(7) 


0.8992(7) 


0.903(1) 


0.748(4) 


43(2) 


0.115(6) 


0.9559(8) 


0.921(5) 


0.880(1) 


0.965(1) 


0.943(2) 


0.806(4) 


58(3) 


0.125(6) 


0.9519(8) 


0.958(3) 


0.877(1) 


0.957(1) 


0.945(2) 


0.817(3) 


61(2) 


0.120(4) 


0.957(1) 


0.949(2) 


0.882(2) 


0.964(2) 


0.949(1) 


0.820(4) 


76.5(6) 


0.129(1) 


0.9532(4) 


0.976(1) 


0.8814(6) 


0.9563(6) 


0.9489(5) 


0.841(2) 


35.0(3) 


0.150(1) 


0.8860(5) 


0.943(2) 


0.8112(7) 


0.8856(6) 


0.8883(7) 


0.761(3) 


43.1(5) 


0.150(2) 


0.8880(2) 


0.945(1) 


0.812(1) 


0.8873(4) 


0.8910(6) 


0.774(5) 


15.0(4) 


0.144(2) 


0.8887(7) 


0.918(7) 


0.813(2) 


0.892(2) 


0.889(1) 


0.69(2) 


12.1(5) 


0.140(1) 


0.8974(3) 


0.920(6) 


0.822(2) 


0.9015(7) 


0.896(1) 


0.70(2) 



Table 7: A chiral extrapolation for the interpolated renormalization factors in the MS-schcmc 
^MS.chir w j t j 1 correS p 0n di n g AWI mass toawi & t momentum transfer /i = 2 GeV. 



m AW i[MeV ] 


a[fm] 


z q 


Z S 


Zy 


Zt 


Z A 


gSub 


31(4) 


0.147(18) 


0.925(1) 


0.969(5) 


0.8453(9) 


0.925(1) 


0.928(2) 


0.799(4) 


43(2) 


0.115(6) 


0.957(1) 


0.891(6) 


0.878(2) 


0.967(2) 


0.940(3) 


0.849(4) 


58(3) 


0.125(6) 


0.962(1) 


0.956(4) 


0.884(2) 


0.967(2) 


0.954(2) 


0.868(3) 


61(2) 


0.120(4) 


0.965(1) 


0.939(3) 


0.886(2) 


0.972(2) 


0.955(2) 


0.867(4) 


76.5(6) 


0.129(1) 
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0.971(2) 
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0.896(2) 


35.0(3) 


0.150(1) 


0.9091(5) 


0.959(2) 


0.8267(7) 
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0.9061(8) 


0.806(2) 


43.1(5) 


0.150(2) 


0.9115(3) 


0.960(1) 


0.8314(8) 


0.9107(4) 


0.9134(7) 


0.819(5) 


15.0(4) 


0.144(2) 


0.9124(6) 


0.941(7) 


0.834(2) 


0.916(1) 


0.912(2) 


0.72(2) 


12.1(5) 


0.140(1) 


0.9202(4) 


0.942(6) 


0.842(1) 


0.9237(7) 


0.919(1) 


0.72(2) 



Table 8: A chiral extrapolation for the reduced, interpolated renormalization factors in the MS- 
scheme jj MS ' chlr w j^ n ^ ne corresponding AWI mass toawi at momentum transfer \i = 2 GeV. 
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Figure 9: The interpolated renormalization factors in the MS-scheme Z MS and Z MS at momentum 
transfer \i = 2 GeV plotted against the lattice spacing a for massive quarks (black) and a chiral 
extrapolation (gray). 
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